Mixing at the external boundary of a submerged turbulent jet 
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We study experimentally and theoretically mixing at the external boundary of a submerged tur- 
bulent jet. In the experimental study we use Particle Image Velocimetry and an Image Processing 
Technique based on the analysis of the intensity of the Mie scattering to determine the spatial 
distribution of tracer particles. An air jet is seeded with the incense smoke particles which are 
characterized by large Schmidt number and small Stokes number. We determine the spatial distri- 
butions of the jet fluid characterized by a high concentration of the particles and of the ambient 
fluid characterized by a low concentration of the tracer particles. In the data analysis we use two 
approaches, whereby one approach is based on the measured phase function for the study of the 
mixed state of two fluids. The other approach is based on the analysis of the two-point second-order 
correlation function of the particle number density fluctuations generated by tangling of the gradient 
of the mean particle number density by the turbulent velocity field. This gradient is formed at the 
external boundary of a submerged turbulent jet. We demonstrate that PDF of the phase function of 
a jet fluid penetrating into an external flow and the two-point second-order correlation function of 
the particle number density do not have universal scaling and cannot be described by a power-law 
function. The theoretical predictions made in this study are in a qualitative agreement with the 
obtained experimental results. 



PACS numbers: 47.27.wg, 47.27.tb 



I. INTRODUCTION 

Jets are often used in numerical, experimental and the- 
oretical studies for modelling of mixing in various engi- 
neering applications (see, e.g., [H, 0, H, 0, H, H, 0, H, H, 
EjIIH)- Different kinds of structures are formed during 
mixing in turbulent jets (see, e.g., [H, H^l)- One of the 
most significant results is a detection of a sharp increase 
of mixing rates caused by the onset of small-scale turbu- 
lence within large-scale coherent motions ^14] . This effect 
that was studied experimentally and numerically, demon- 
strates complex nonlinear dynamics of mixing (see, e.g., 
BIO, EH [13, [Hi)- However, it is not clear how such 
mixing states are attained, and the ambiguity is exacer- 
bated by the differences in experimental results. 

The possible reason for such variety of experimental re- 
sults is the difference in Schmidt number between gases 
and liquids. It was found in 17J| that mixing becomes 
more effective in a liquid jet for the larger Reynolds 
numbers. However, the dependence of the normalized 
scalar variance on Reynolds number for the gaseous jet 
is very weak. It was suggested in ^7|] that this diffcr- 
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ence is caused by the Schmidt number effect. In partic- 
ular, when Schmidt number Sc ~ 1 the larger molecular 
species diffusivity results in better mixing of the scalar 
field even at lower Reynolds numbers. Here Sc= v/Dm 
is the Schmidt number. Dm is the coefficient of Brownian 
(molecular) diffusion and v is the kinematic viscosity. For 
the liquid jets, Sc ~ 10^, the improved mixing requires 
the enhancement in the interfacial surface-to- volume ra- 
tio (smaller distances between iso-scalar surfaces) asso- 
ciated with higher Reynolds numbers. Due to the large 
differences in Schmidt numbers the results of these stud- 
ies in liquids are not directly applicable to the case of 
gaseous shear-layer mixing. 

In the present study we investigate experimentally 
mixing at the air jet interface using the incense smoke 
that is characterized by significantly larger Schmidt num- 
ber (Sc~ 10^) than that employed in the previous stud- 
ies of gas flow mixing. The high magnitude of Schmidt 
number prevents from fast molecular diffusion to affect 
mixing. We study mixing at the external boundary of a 
submerged turbulent jet using Particle Image Velocime- 
try (PIV) and an Image Processing Technique based on 
the analysis of the intensity of the Mie scattering in order 
to determine the spatial distribution of the jet fluid (with 
a high concentration of the tracer particles) and of the 
ambient fluid (with a low concentration of the tracer par- 
ticles) . In this flow there is no a sharp interface between 
the jet and ambient air. 

In the data analysis we use two approaches, whereby 
one approach is similar to that used previously in the 
analysis of Ray leigh- Taylor instability (see [l^). This 
approach is based on the measured phase function for 
the study of the mixed state of two fluids. The other ap- 



2 



proach used in our study is based on the analysis of the 
two-point second-order correlation function of the parti- 
cle number density fluctuations generated by tangling of 
the gradient of the mean particle number density by the 
turbulent velocity field (see, e.g., 0, HH). This gradi- 
ent is formed at the external boundary of a submerged 
turbulent jet. Both approaches demonstrate that proba- 
bility density function (PDF) of the phase function of a 
jet fluid penetrating into an external flow and the two- 
point second-order correlation function of the particle 
number density do not have universal scaling and cannot 
be fitted by a power-law function. There is a qualitative 
agreement between theoretical and experimental results 
obtained in our study. 

The paper is organized as follows. In Sect. II we de- 
scribe the experimental set-up for a laboratory study of 
mixing at the external boundary of a submerged tur- 
bulent jet. In Sect. Ill we discuss the experimental re- 
sults and their detailed analysis by means of the approach 
based on the measured phase function. In Sect. IV we 
perform a theoretical study of mixing at the external 
boundary of a submerged turbulent jet based on the anal- 
ysis of the two-point second-order correlation function of 
the particle number density. In Sect. IV we also com- 
pare the theoretical predictions with the obtained exper- 
imental results for the correlation function of the par- 
ticle number density. Finally, conclusions are drawn in 
Sect. V. In Appendix we perform an asymptotic analysis 
of the solution of equation for the two-point second-order 
correlation function. 



II. EXPERIMENTAL SET-UP 

The experimental set-up with air jet shown in Fig. [1] 
includes the chamber with the transparent plexiglass 
walls (1), and the cylindrical tube (2) with a jet nozzle 
mounted at the tip. A submerged air jet (3) discharges 
into the chamber. The cross-section of the channel with 
transparent walls is 47 x 47 cm^, a diameter of the jet 
nozzle is = 10 mm and the diameter of the jet in 
the probed region is about 35 mm. Consequently, the 
side walls of the channel weakly affect the development 
of the jet. The optics (4) produces a fight sheet (5), 
and the scattered light from the probed flow region (6) is 
recorded with the CCD camera (7), see Fig.[TJ The laser 
sheet thickness is 1 mm. Note that in the present paper 
we study mixing of particles at the external boundary of 
a jet and we do not investigate the dynamics of jet. 

Flow measurements have been conducted with digital 
Particle Image Velocimetry (PIV) using a LaVision Flow 
Master III system to determine the jet velocity field (see, 
e-g-, QlL i23„ 2M)- The fight sheet is provided by a dou- 
ble PIV Nd-YAG laser. Continuum Surelite 2 x 170 mJ. 
The beam of the laser is directed with the aid of a few 
prisms to the light sheet optics that comprises spherical 
and cylindrical Galilei telescopes with tunable divergence 
and adjustable focus length. A progressive-scan 12 Bit 




FIG. 1; A scheme of the test section and PIV technique; a 
sketch of a jet coordinates and a range of measurements. 1 - 
channel with transparent walls, 2 - tube with a jet nozzle, 3 - 
submerged jet, 4 - light sheet optics, 5 - laser light sheet, 6 - 
image area, 7 - CCD camera. 

digital CCD camera with dual-frame-technique for cross 
correlation captures the images of 1280 x 1024 pixels with 
a size 6.7 /im x6.7/im. A programmable Timing Unit 
(an interface card of the system PC) generates sequences 
of pulses to control the laser, the camera and data ac- 
quisition. The software package DaVis 7 is applied to 
control all hardware components and for image acquisi- 
tion and visualization. Jet flow images are analyzed with 
the package supplemented with the software developed 
in this study. 

In this experimental study we use methods of measure- 
ments applied previously in different studies of turbulent 
transport of particles (for details see [H, [H, [131). In 
particular, we use PIV and an Image Processing Tech- 
nique based on the analysis of the intensity of the Mie 
scattering to determine the spatial distribution of tracer 
particles (see, e.g., [III). As a tracer we use smoke pro- 
duced by incense particles sublimation. The submerged 
air jet, seeded with the incense smoke particles of sub- 
micron sizes with the mean diameter 0.7/^m, flows into 
the still air inside the chamber. An estimate of the mean 
response time of these particles Tp to a step change of the 
velocity is of the order of 10~^ s. Since the latter value is 
much smaller than the Kolmogorov time of the flow 
in the jet (r^ varies in the rage from 10~^ s to 10"'^ s), 
and smoke particles Stokes number St= Tp/T^^ <C 1, it 
can be safely assumed that these particle follow the flow 
faithfully. 

Laser sheet light scattered by the incense smoke parti- 
cles is captured with a CCD camera. Series of 50 pairs of 
images with the adjustable time delay are acquired with 
a frequency of 2 Hz and are stored for ensemble and spa- 
tial averaging of flow characteristics. We measure spatial 
particle distribution in a flow area of 18.4 x 18.4 cm^ 
with a spatial resolution of 1024 x 1024 pixels that is 
0.18 mm/pixel. 

In the experiments the jet fluid is characterized by a 
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high concentration of the tracer particles, and the ex- 
ternal fluid is characterized by a low concentration of 
the particles. Every recorded image is normalized by 
a light intensity measured just at the jet entrance into 
the chamber in order to eliminate the effects associated 
with a change of concentration of the incense smoke. We 
measure parameters of mixing in a local system of coor- 
dinates (y, z), where y-axis is directed along the averaged 
jet-ambient fluid interface, z-axis is perpendicular to this 
interface (see Fig.[T]). The jet parameters are determined 
in intervals L, where the flow characteristics are statisti- 
cally homogeneous. 



III. DATA PROCESSING BASED ON THE 
MEASURED PHASE FUNCTION 

In this section we use the approach applied previously 
in the analysis of Ray leigh- Taylor instability (see p^). 
This approach is based on the measured phase function, 
X{y,(j)), defined as follows: X{y,(j)) — 1 when the point 
{y, (f) is within the jet fluid and X{y, (f)) — Q elsewhere 
(see [2^), where (j) is the polar angle (see Fig. [1]). We 
analyze the statistical properties of the system by aver- 
aging over an interval of L along the y direction at fixed 
angles 0, i.e., {X) = J^^ X {y , (j)) dy , where L is the 
length of the interval in y direction (see Fig. [T]). 

Let us introduce the first and the second moments of 
the measured phase function: «(</>) = (X) and ?(A, cj)) = 
{X{y - A/2,(/))X(y + A/2,0)) , where a((/)) character- 
izes the fraction of y which is occupied by the jet fluid, 
A are the sizes of regions (with X{y) = 1) occupied by 
the jet fluid (see Fig. (2). The study of the mixed state 
of two fluids induced by Rayleigh- Taylor instability con- 
ducted in [19] indicates that PDF of the random size A 
of the regions occupied by the jet fluid can be described 
by Gamma distribution 



/(A) 



1 aL exp(-A/A) 
r(2-r) A2" {A/xy 



(1) 



where A are the sizes of the regions occupied by the jet 
fluid, the scale A and the exponent r are parameters char- 
acterizing the length scale and a deviation from the ex- 
ponential PDF. A ratio of the moments 
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determines a characteristic scale A„ , where n is the num- 
ber of the statistical moment. Equation ([1]) for PDF im- 
plies that 



— = const . 
n 



(3) 



This property of the PDF of the sizes A of the regions 
occupied by the jet fluid, is valid if the parameter A in 
Eqs. HI) and ^ varies in the infinite interval < A < cx). 
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FIG. 2: The measured phase function, X{y) at fixed angle 
(j). Here A are the sizes of regions with X{y) = 1 that are 
occupied by the jet fluid. 



We show below that the measured PDF of the phase 
function in our experiments can be reasonably approx- 
imated by Gamma distribution. In the present analy- 
sis we use images of a jet obtained in our experiments, 
whereby the interface between a jet and an ambient fluid 
is determined for a mean jet image averaged over 50 in- 
dependent instantaneous images. The images are con- 
verted into a binary form using the following procedure. 
The threshold for the image transformation into a binary 
form is chosen equal to the intensity whereby the PDF 
of the scattered light intensity in the jet equals to that in 
the ambient fluid (see Fig. [H where this threshold is 96 
counts). Thereafter all pixels of every image where the 
scattered light intensity is above the threshold value, are 
assigned the value of the phase function X(?/, 0) equal 
1 and are categorized as a jet fluid, while pixels with 
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FIG. 3: Measured PDF of light intensities in the jet (right 
hne) and in the ambient fluid (left line). 



light intensity below the threshold are assigned the value 
of the phase function equal and are categorized as an 
ambient fluid. 

After implementing this procedure the mean image of 
the instantaneous binary jet images is obtained. The 
mean image is transformed into a binary form with the 
threshold 0.5. In particular, the pixels, where a proba- 
bility of a jet fluid appearance was higher then 0.5 were 
assigned the value of the phase function equal 1, while 
all the rest were assigned the value of the phase func- 
tion equal 0. The boundary of a jet is a line, where a 
mean value of a phase function equals 0.5. The threshold 
adopted in this study (the intensity whereby the PDFs 
of the scattered light intensity in the jet equals to that 
in the ambient fluid), allows to decrease a systematic er- 
ror during the transformation into a binary form. Our 
analysis has shown that the obtained results on mixing 
of particles at the external boundary of a jet are not very 
sensitive to a particular choice of the threshold. 

We can roughly distinguish between three regions in 
the jet down steam the nozzle. In the first region having 
the size ~ D the jet is stable, where D is a diameter of 
a nozzle. The second region includes an unstable part of 
the jet, whereby the Kelvin-Helmholtz instability might 
be the main feature of the flow, and this region extends 
up to a distance of about ~ 313. The third region down 
steam is a turbulent flow region that is the subject of our 
study. 

A phase function, X{y, z), is determined for each image 
in the interval oly/D from 4.5 to 10.9 along the different 
directions with a given (j) (or for different values of z in the 
range ±Z3), see Fig.[T] We determine a mean value, a, of 
the phase function X(y, z) averaged over the ensemble of 
the jet images. For instance, the function a{z) is shown in 
Fig. [4] for two experiments with different global Reynolds 
numbers, 7^ = 8.4 x 10^ and 7^ = lO-*. Here n = VD/v 
is the global Reynolds number based on the mean jet 
velocity V at the nozzle exit and the diameter Z) of a 




FIG. 4: First moment a of the phase function for different z 
inside and outside the jet determined in the experiments with 
different global Reynolds numbers: TZ = 8.4 x 10'^ (triangles) 
and TL = 10* (squares). Here 2 = corresponds to the mean 
boundary of the jet, z < inside the jet and z > outside 
the jet. 




-1.4 -1.2 -1 -0.8 -0.6 -0.4 logj|^(A/D) 

FIG. 5: PDF of the phase function of the jet fluid pene- 
trating into an ambient flow (squares) and its approximation 
(solid line) for z = by the Gamma distribution (see Eq. ^) 
with the following parameters: r — 1.03, \/D = 0.13, and 
A — 0.182. For comparison, a power- law fitting curve Ai 
(dashed line) is also shown, where ri = 2.46, X/D = 0.13, and 
Ai = 3.36 X 10^^ 

jet nozzle. Although we performed experiments for two 
values of the Reynolds number, we did not study the 
dependence of mixing on the Reynolds number. The line 
z = corresponds to the mean boundary of the jet, z < 
inside the jet and z > outside the jet. The mean 
value a of the phase function decreases along the positive 
direction of the axis z. 

In our analysis of mixing of particles at the external 
boundary of a jet it is not important whether the jet 
at distances y/D — 4.5 — 10.9 could be affected by the 
initial conditions. Indeed, the interface between the jet 
and the ambient fluid is determined for a mean jet image 
averaged over 50 independent images that are acquired 
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FIG. 6: Ratios of the moments \n/nD of the measured PDF 
versus the length A^/f for n varied from 1 (squares) to 4 (cir- 
cles), where Am is the maximum size of the observed jet fluid 
inclusions. Ratios A„/n_D are also determined using a direct 
fit of the measured PDF to Gamma distribution [Eqs. ((T}- 
([2])] with the following parameters: r — 1.03, X/D — 0.13, 
and A = 0.182, where n = 1 is the upper solid curve, and 
n = 4 is the lower solid curve) . 



with a frequency of 2 Hz. The successive pairs of images 
are not correlated since the turbulence correlation time 
To is of the order of 10~^ s and the characteristic time 
of the mean flow in the image scale is of the order of 
10~^ s. The estimate for the turbulence correlation time 
is obtained using the measured characteristic turbulent 
velocity mq ^ (1 — 2) x 10"^ cm/s, the maximum scale of 
turbulent motions ^ 1 cm and tq = Iq/uq ~ 10"'^ s. 

PDF of the phase function of a jet fluid penetrating 
into an ambient fluid is determined for the same intervals 
oiy/D for different z using parametric histogram estima- 
tor (see, e.g., [30[). The obtained PDF is approximated 
by the Gamma distribution ([T]) using three adjustable 
parameters r, A and A — aL/T{2 — r). The measured 
in our experiments PDF of the phase function and its 
approximation by the Gamma distribution, are shown in 
Fig. [51 where the parameters of the Gamma distribution 
are r = 1.03, \/D = 0.13, and A = 0.182. These PDFs 
vary more slowly for large values of A than for small 
A. For comparison, a power-law fitting curve Ai ^^^i is 
also shown in Fig. [5l where ri = 2.46, X/D — 0.13, and 
Ai = 3.36 X 10~^. Clearly, the measured PDF cannot be 
fitted by a power-law function (see Fig. [5]). 

The ratios of the moments A„/ri for the measured PDF 
versus the length A.,„ are shown in Fig. [HI where Am is 
the maximum size of the observed jet fluid inclusions. A 
range of our experimental data is limited to Am/D = 0.6, 
where the probability of the event with A/D = 0.6 is 
about 3 X 10^"'. Consequently, it can be expected that 
only one such event occurs in a sample with the size 
of about 3300 events (not to mention the events with 
A/D > 0.6). Therefore, it is not feasible to expand sig- 
nificantly the Am/ D range of the observed events, since 



FIG. 7: Dependence of the exponent r on the distance from 
the jet boundary for 7^ = 8.4 x 10^ (triangles) and 7^ = 10* 
(circles). 
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FIG. 8: Dependence of a scale A on a distance z/D from 
the jet boundary for 7^ = 8.4 x 10^ (triangles) and 7^ = 10* 
(circles). 



a probability of the events with A/D > 0.6 is extremely 
low. Using Eq. ([2]) we also determine the scale A„/n 
for n varied from 1 to 4 using a direct fit (solid lines in 
Fig. [5]) of the measured PDF to Gamma distribution with 
the following parameters: r = 1.03, X/D — 0.13, and 
A — 0.182. Inspection of Fig. [6] shows that if the range 
of integration Am is too short (i.e., A,„ < 10A„/n), the 
ratios Xn/n are different for different n. On the other 
hand, when Am > 10A„/n, the ratios Xn/n are inde- 
pendent of n. Therefore, the method of determining the 
parameters of PDF [see Eq. (IT|)], using the property ([3]) 
of Gamma distribution, A„/ri=const, is not appropriate 
if A,„ < IQXn/n. The latter condition corresponds to 
all our experiments (see Fig. [6]) . This implies that in our 
experiments we cannot account for large A when we de- 
termine A„. Therefore, the method based on Eq. ([3]) can 
be used only when Am > 10 Xn/n. 
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The parameters r and A of the PDF of the phase func- 
tion versus z determined in the experiments are shown in 
Figs. [7] and [H The characteristic scale A that is shown in 
a logarithmic scale (see Fig. varies significantly, de- 
creasing by 1.5 orders of magnitude from the internal to 
the external regions of a jet. This dependence is close to 
an exponential function. The behaviour of the functions 
\{z) and a{z) is similar (see Figs. S] and [51). Indeed, a 
is the first moment (mean value) of the phase function, 
while A should coincide with the first moment Ai if the 
sizes A of the regions occupied by the jet fluid vary from 
to very large values. For negative z (i.e., inside the 
jet) the size A of the region occupied by the jet fluid is 
much larger than that for positive z (i.e., outside the jet), 
whereby the concentration of the jet fluid is small. 

It must be noted that we perform experiments in the 
region located at the distance 4.5 — 11 cm downstream 
from the nozzle. The reason is that in this region the gra- 
dient of the mean particle number density at the external 
boundary of a submerged turbulent jet is very large. In 
this region the tangling mechanism that is responsible 
for particle mixing, is more pronounced. The tangling 
mechanism depends on the gradient of the mean particle 
number density (see Sect. IV). Increase of the distance 
from the nozzle, is accompanied by a decrease of the gra- 
dient of the mean particle number density. In addition, 
in the regions located far from the nozzle, the effect of 
the bottom wall of the chamber (e.g., appearance of the 
secondary flows) becomes stronger. 

Note that mixing in water jet flows was also studied in 
recent experiments [3l| (see also [H,!!^). In these exper- 
iments the water jet with ink having a diameter D = 0.8 
cm and the global Reynolds numbers TZ ^ 10^, discharges 
into a square duct with a width L = i cm. The mea- 
surements described in [sH were performed in the duct 
section at the distance x/D in the range from 6.6 to 22.5 
downstream of the jet entrance (see Figure 2b in [sH]). 
Inspection of Figure 1 in [3ll | reveals that the mixed part 
of the flow spreads over the whole cross-section of the 
channel after ~ 4.5 cm {x/D ~ 5.6) downstream of the 
jet entrance. Consequently, most of the measurements in 
[3l| were performed in the non-turbulent duct flow with 
the global Reynolds number 7?. ~ 2 x 10"^ (based on the 
width of the duct and the mean flow velocity). The flow 
in another experiment reported in [3lj | is also laminar 
{TZ ^ 0.1). In these experiments it was found that the 
PDF of the concentration of the dye in the channel can 
be described by Gamma distribution with a positive pa- 
rameter r = n — 1~3 — 84 in a power function. In our 
study we consider a kinematic problem of mixing at the 
turbulent fluid-ambient fluid interface and found that the 
PDF of the random size of the regions occupied by the 
jet fluid can be described by a Gamma function with a 
negative parameter r ~ —0.9 in a power function. These 
flndings imply the absence of power-law scaling. 

Although the PDF of phase function was fltted in Fig.[S] 
by Gamma function, one of the main flndings of this sec- 
tion is that the PDF of phase function does not have a 



universal scaling and cannot be described by a power-law 
function. 



IV. THEORETICAL MODELLING OF 
TURBULENT MIXING AND COMPARISON 
WITH EXPERIMENTAL RESULTS 

In the previous section we analyze experimental data 
using the measured phase function. In this section we use 
another approach in the analysis of turbulent mixing of 
particles at the external boundary of the submerged tur- 
bulent jet. This approach is based on the analysis of the 
two-point second-order correlation function of the parti- 
cle number density fluctuations generated by tangling of 
the gradient of the mean particle number density by the 
turbulent velocity fleld. This gradient is formed at the 
external boundary of a submerged turbulent jet. 

The equation for the number density nit, x) of particles 
advected by a fluid velocity fleld reads: 



dn 

— + W-{nv) ^ D„,An 
ot 



(4) 



where Dm is the coefficient of molecular (Brownian) dif- 
fusion, v(t, x) = V -I- u is the particle velocity field, 
V = (v) is the mean velocity and u are velocity fluc- 
tuations. Equation ^ implies conservation of the total 
number of particles in a closed volume. In this study we 
consider incompressible particle turbulent velocity fleld 
(V-u = 0). This implies that we neglect a small com- 
pressibility of particle velocity fleld caused by particle 
inertia, because the particle size in our experiments is 
less that 1/im. We also neglect very small compressibil- 
ity of the air velocity. Consider 0(t,x) ^ n — N, the 
deviation of n from the mean number density of particles 
N. The equation for fluctuations of the particle number 
density 6(t,x) reads: 

^ + v-(eu-(eu)) = i?™Ae-(u.v)iv, (5) 

(see, e.g., [S^, HI])- Using Eq. ^ we derive equation 
for the evolution of the two-point second-order correla- 
tion function of the particle number density, $(t,R) = 
(e(t,x)e(t,x + R)): 



'dt 



^ 2 



A,(o)-A,(R) 



V^Vo$-l-/ 



(6) 



where Dy(R) is the turbulent diffusion ten- 
sor (see below), Sij is the Kronecker tensor, 
/ = (ViV)^ exp(— c*i?/Zo) is the source of par- 
ticle number density fluctuations at the jet-ambient 
fluid interface, D^ oc Iq uq is the turbulent diffusion 
coefficient, uq is the characteristic turbulent velocity in 
the maximum scale Iq of turbulent motions, > 1 is a 
free constant and the angular brackets denote ensemble 
averaging. The source function / is related to the last 
term in the right hand side of Eq. 



7 



When VA'' ^ the nonzero source / results in gen- 
eration of fluctuations of the particle number density 
caused by tangling of the gradient of the mean particle 
number density by the turbulent velocity field. In this 
section we study this effect in detail. Equation ([6]) has 
been derived in [sH] for a delta-correlated in time ran- 
dom velocity field. For a turbulent velocity field with a 
finite correlation time Eq. ^ has been derived in [2l| 
by means of stochastic calculus, i.e., Wiener path inte- 
gral representation of the solution of the Cauchy problem 
for Eq. ([?]), using Feynman-Kac formula and Camcron- 
Martin-Girsanov theorem. The comprehensive descrip- 
tion of this approach can be found in [13, HU, [131 . 

T 

The turbulent diffusion tensor, D,y(R), is given by the 
following formula: 



/•OC 

dI = / {u^ [0, ^{t, x|0)] u, [r, ^{t, X + R|t)] ) dr 
Jo 



(7) 



(sec "21^). Here the Wiener trajectory ^(<,x|s) (which is 
often called the Wiener path) is defined as follows: 



|(i,x|s) = X - / u[T,^{t,^\T)] dT~ ^2D^Vf{t 
J s 



(8) 



where w(i) is the Wiener random process which describes 
the Brownian motion (molecular diffusion). The Wiener 
random process w(i) is defined by the following proper- 
ties: {■w{t))^^0, {wi(t + T)wj(t))^ = rSij, and (. . .)w 
denotes the mathematical expectation over the statistics 
of the Wiener process. The velocity Mi[T, ^(t, x|r)] de- 
scribes the Eulerian velocity calculated at the Wiener tra- 
jectory. For a random incompressible velocity field with 
a finite correlation time the tensor of turbulent diffusion 
isgiven by dI{R) w T-\^,{t,x\0) ^j{t,x + R\t)) (see 
|2l|). Hereafter (...) denotes averaging over the statistics 
of both, turbulent velocity field and the Wiener process. 

If the turbulent velocity field is not delta-correlated 
in time (e.g., the correlation time is small yet finite), 
the tensor of turbulent diffusion, £)jj(R), is compressible, 
i.e., {d/dR,)D^-{R) ^ (for details see [2l|). Let us 
consider the parameter tr^ that characterizes the degree 
of compressibility of the tensor of turbulent diffusion: 



V,V,Z?,,(R) 



((VX|)2) 



(9) 



where e^-fc is the fully antisymmetric Levi-Civita unit ten- 
sor. Note that when the turbulent velocity field is a delta- 
correlated in time random process, V-^ oc V-u, while for 
a turbulent flow with the finite correlation time V-^ has 
a contribution that is independent of V-u. 

Let us discuss the assumptions underlying the above 
model of particle transport in turbulent flow. We use 
the tensor of turbulent diffusion D^^ (R) for isotropic and 



homogeneous turbulent flow. In our experiments the ve- 
locity field is generally anisotropic. We determine the 
two-point second-order correlation function <I>(t, R) in 
the vicinity of the mean boundary of the jet along y- 
axis, so that the size of the probed region in z direction 
is small. It is of the order of the maximum scale of tur- 
bulent motions Zo- 

In addition, a main contribution to the level of fluc- 
tuations of particle number density is due to the mode 
with the minimum damping rate (see Eqs. ([22]) and ((23|) 
below). This mode is an isotropic solution of Eq. ([6|. 
Consequently, it is plausible to neglect the anisotropic ef- 
fects. This assumption is also supported by our measure- 
ments of the two-point second-order correlation function 
$(t, R) determined along y-axis for different values of z 
inside and outside the jet (see Fig. [T51 below). The differ- 
ence between these correlation functions determined for 
different z is very small. We perform measurements in 
the range along y-axis whereby the turbulence is nearly 
homogeneous. On the other hand, the turbulence param- 
eters Zq and Wo vary slowly with z in the probed region. 

Note that the mechanism of mixing related to the tan- 
gling of the gradient of the mean particle number den- 
sity at the external boundary of a submerged turbulent 
jet by the turbulent velocity field, is enough robust. The 
properties of the tangling are not very sensitive to the 
exponent of the energy spectrum of the background tur- 
bulence. The requirements that turbulence should be 
an isotropic, homogeneous, and should have a very long 
inertial range (a fully developed turbulence) , are not nec- 
essary for the tangling mechanism. 

Anisotropy effects can complicate the theoretical anal- 
ysis, but do not introduce new physics in the mixing pro- 
cess. The reason is that the main contribution to the 
the tangling mechanism is at Kolmogorov (viscous) scale 
of turbulent motions. At this scale turbulence can be 
considered as nearly isotropic, while anisotropy effects 
can be essential in the vicinity of the maximum scale of 
turbulent motions. 

T 

Using these arguments, we consider the tensor D^^ (R) 
in the following form: 



D^AR) = [F{R) + F,{R)]5,j + RF^ 



I RiRj 



i?2 



RF' 



RiRj 
~R2- 



(10) 



IoUq/3 is the turbulent diffusion coefficient, 
F(0) 1 - Fc(0) and F' = dF/dR. The function 
Fc{R) describes the compressible (potential) component, 
whereas F(R) corresponds to vortical (incompressible) 
part of the turbulent diffusion tensor. Using Eqs. ^ 
and (|10[) we derive equation for the two-point second- 
order correlation function ^(R) written in a dimension- 
less form: 



9$ 
'dt 



1 



M{R) 



$" + 2(i+x(i?)) 



(11) 
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where time t is measured in units of tq = Iq/uq, distance 
R is measured in units of Iq, and 



1 



M{R) 

x{R) 



2 2. 
M{R) 



1 - F - {RF,)'] 
{F - 2F,y , 



(12) 
(13) 



and Pe = uqIq/Dj^ is the Peclet number. The two-point 
correlation function $(i?) satisfies the following bound- 
ary conditions: ^{R = 0) = 1 and $(i? ^ oo) = 0. This 
function has a global maximum at i? = and therefore 
it satisfies the conditions: 

$'(i? = 0) = 0, = 0) < , 

= 0) > > 0)1 . 

Let us introduce a function 



V'(<,R) = i?$(i,R)exp 



X{x) dx 



(14) 



Equations (fTTjl and ([T4|) yield the following equation for 
the function ijj{R): 



d'4' 
'dt 



where 



U{R) 



1 d^'ijj 
M{R) 

x{R) 



- U{R)^ + /o 



M{R) 
lo = R I{R) cxp 



x{R) + - 

X 



x{x) dx 



(15) 

(16) 
(17) 



In the analysis we use a quantum mechanics anal- 
ogy, whereby Eq. (fT5|) is regarded as a one-dimensional 
Schrodinger equation with a variable mass. We seek for 
the solution of Eq. (fTS]) in the following form: 



V'(t,i?)=5^0p(t)*p(i?) 



(18) 



where ^p(i?) are the eigenfunctions determined by the 
following equation: 



1 d^^p 
M{R) dR^ 



[27p + [/(i?)]*p-0 



(19) 



The condition of the orthogonality for the eigenfunctions 
reads: 



M(i?)*i(i?)«'j (i?) dR 



(20) 



Substituting solution into Eq. (jlSp. multiplying the 
obtained equation by (i?) and integrating over R yields 
the following equation for the function 4>p{t): 



d 

di^P 



(21) 



In derivation of Eq. (|2T|) we use Eq. ([20|) and take into 
account that 



oc / t'R 

M{R)RI{R) e^pi^J x{x)dx]'^p{R)dR 

2 '2/T7AT\2 



3c, 



where the function AI {R)I{R) is strongly localized and is 
approximated by the delta function. We also assume here 
that the basis of the eigenfunctions ^'p(i?) is complete. 
The solution of Eq. ((2T|l for the function 4>p{t) with the 
initial condition (j)p(t = 0) — reads: 



= ^A(ViV)' [l - exp(-|7p|i) 
oc, |7p| L 



(22) 



Equations (O, (|18p and ([22)1 yield the correlation func- 
tion $(t,i?): 



Ht,R) 



3c^ 



l-exp(-|7p|i) 



(23) 



where we used an estimate (VTV)^ ^ {SNy/Ll. The 
main contribution to the correlation function $(i, R) for 
t:S> |7p|^^ in Eq. (^5]) is due to the mode with the mini- 
mum damping rate |7p|, i.e., for the mode with p = 1. 

Particular formulas for the potential U{R) and the 
mass M{R) depend on the functions F{R) and Fc{R). 
For instance, we may choose these functions in the fol- 
lowing form 



FiR) 
FciR) 

fiR) 



l + (j. 



exp[-/(i?)] , 



— ^ expha,/(i?)] 

1 + (Ty 

R^ 

Re-i/2 + ^2/3 ' 



(24) 
(25) 

(26) 



where Re — Uf)l^/v is the Reynolds number based on 
turbulent velocity uq and maximum scale lo of turbulent 
motion. Equation ()26p is similar to an interpolation for- 
mula derived by Batchelor for the correlation function of 
the velocity field that is valid for a turbulence with Kol- 
mogorov spectrum in the inertial range and for random 
motions in the viscous range of scales (see, e.g., [38l,[39|). 
In particular, in the inertial range of turbulent scales, 
Re"^''^ <C i? <C 1, the correlation function for the tur- 
bulent diffusion tensor is F{R) cx 1 — R^^^, where R is 
measured in units of maximum scale of turbulent mo- 
tions Zq. The corresponding correlation function for the 
turbulent velocity field cx 1 — R^/^ . The difference be- 
tween the scalings of the turbulent diffusion tensor and 
the correlation function of the turbulent velocity field is 
caused by the scaling of correlation time r(i?) oc R^/^. 
On the other hand, in the viscous range, R <C Re"^^'*, 
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1 2 3 R/ljj 

FIG. 9: Potential U versus R/lo for = 1 and different 
values of parameter Uc. ac — 1.2 (dashed-dotted line), Uc — 2 
(thick solid line) and ac — i (dashed line). For comparison 
the case a^, — (the delta-correlated in time velocity field) is 
shown (thin solid line). Here Re= 2 x 10^ and Pe= 2 x 10*. 
Horizontal dashed-dotted line U = 0.09 corresponds to the 
damping rate for the first mode 71 = — 0.09r(J~^. 



the correlation function for the turbulent diffusion tensor 
is F{R) cx 1 — Re^^^i?^. This correlation function is sim- 
ilar to that for the velocity field because in the viscous 
range the correlation time is independent of scale. On 
the other hand, for large scales I there is no turbu- 
lence, so that for R > 1 the functions F{R) and Fc{R) 
should decrease sharply to zero. 

The particular choice of the functions F{R), Fc{R) 
and f{R) in the paper describes the well-known proper- 
ties of turbulent velocity field obtained from laboratory 
experiments, numerical simulations and theoretical stud- 
ies (see, e.g., [sl, HI]). In this study we choose the ex- 
ponential form for the functions F{R) and Fc{R). The 
final results are not sensitive to the form for these func- 
tions at large scales i? ^ 1 . Note also that the particular 
choice of the function f{R) is also not very important. 
It should include turbulent motions in the inertial range 
(e.g., a turbulence with Kolmogorov spectrum in the in- 
ertial range) and random motions in the viscous range of 
scales. 

Using Eqs. we determine the potential U{R) 

and the mass M{R) in Eq. (fTQ]). For instance, the poten- 
tial U (R) for different values of parameter Oc is shown 
in Fig. [HI This potential U{R) for smaller scales is also 
shown in Fig. [TOl The function 1/M{R) for different 
values of parameter ac is shown in Fig. 1111 Very impor- 
tant features of the potential U (R) which determine the 
asymptotic behaviour of the solution of Eq. p^ . are as 
follows. At small scales, the potential U{R) is negative 
for ^ 0, and there is a positive maximum of the po- 
tential U{R) at the scales i? ~ 1. On the other hand, 
when the degree of compressibility of the turbulent diffu- 
sion tensor cr^, = (e.g., for the delta-correlated in time 




.01 .1 1 10 100 R/ri 

FIG. 10: Potential U/Ke^^^ versus R/rj at small scales for 
(7y = 1 and different values of parameter Uc- = 1.2 

(dashed-dotted line), ac — 2 (thick solid line) and ac = 4 
(dashed line). For comparison the case cr^ = (the delta- 
correlated in time velocity field) is shown (thin solid line). 
Here tj = lo/Re'^^* is the Kolmogorov (viscous) scale and 
Tr] = ro/Re^''^ is the Kolmogorov time, Re= 2 x 10'^ and 
Pe= 2 X 10^ 



1/M(R) 




FIG. 11: 1/M versus R/lo for a^, — 1 and different values of 
parameter ac'. ac — 1.2 (dashed-dotted line), ttc = 2 (thick 
solid line) and ttc = 4 (dashed line). For comparison the 
case CTy = (the delta-correlated in time velocity field) is 
shown (thin solid line). The horizontal thin dashed-dotted 
line determined by equation, = 2/3, corresponds to the 

scales at which there is no turbulence. Here Re= 2 x 10"^ and 
Pe= 2 X 10*. 



turbulent velocity field), the potential U{R) is positive 
for all scales and decreases for the larger scales, so that 
there is no maximum of the potential U (R) for R > 0. 
Analysis performed by Kraichnan in j36j for the delta- 
correlated in time turbulent velocity field, showed that 
the characteristic damping rate of the particle number 
density fluctuations is very high 7 ~ — ^/Rstq^, 

where = To/V^e is the Kolmogorov time (i.e., the 
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turbulent time at the viscous Kolmogorov scale). The 
latter implies that the level of these fluctuations is very 
low - (7ro)-i - Re"i/2 (ggg 

In a real flow with a finite correlation time the degree 
of compressibility of the turbulent diffusion tensor ct^ ^ 
0, the characteristic damping rate of the fluctuations of 
the particle number density is not high 7 < Tq^, and 
therefore, the level of these fluctuations is not small. In 
this case the maximum of the potential U = J7,nax is 
attained at the scale i?rnax 1, whereby the inverse mass 
^ 1 and changes with R considerably slower than 
that in the range i? <C 1. The asymptotic solution for 
the two-point second-order correlation function in 
this region has a form of the Gamma distribution ^{R) cx 

exp[-K(i?- i^max)]/-R^ whcrc K = (C/inax - 2|7|)^^^ and 
6 ~ 1. The detailed asymptotic analysis of the solution 
for the two-point second-order correlation function ^{R) 
in different ranges of scales is performed in Appendix. 

Now we solve Eq. (jl9p numerically in order to de- 
termine the two-point second-order correlation function 
^p{R). The normalized correlation function for the first 
mode (p = 1) for cr^ = 1 and different values of param- 
eter ttc is shown in Fig. [121 Note that the parameter Oc 
strongly affects the compressible part of the tensor of tur- 
bulent diffusion caused by the finite correlation time of 
turbulent velocity field. The location of the maximum in 
the potential U{R) strongly depends on the parameter Oc 
and is nearly independent of the Reynolds number. The 
damping rate for the first mode 71 = —0.09tq^. The 
numerical solution for the correlation function <l>(i?) is 
in agreement with the results of the asymptotic analysis 
performed in Appendix. 

In order to compare the theoretical predictions with 
the experimental results, the analysis of the experimental 
data is performed also without transformation of images 
into a binary form. In particular, the two-point second- 
order correlation function ^{R) determined in our exper- 
iments without and with transformation of images into a 
binary form is shown in Fig. 1121 The distance between 
two points R in our experiments is measured along the y 
axis at 2; = (at the mean boundary of the jet) . In Fig. [13] 
we also show the correlation functions determined in our 
experiments (without transformation of images into a bi- 
nary form) for different z inside and outside the jet. 

In Figs. [T^ and [T51 we use the Reynolds number. Re, 
based on the turbulent velocity velocity uq and the tur- 
bulent correlation length ^o, where the characteristic tur- 
bulent velocity uq ^ Q.2V and ^ D, so that Re= 0.27?.. 
Note that the global Reynolds number TZ based on the 
mean jet velocity V at the nozzle exit and the diameter D 
of a jet nozzle, is more often used in experimental studies 
(e.g., for comparisons of different experimental set-ups), 
while the Reynolds number. Re, based on the turbulent 
velocity velocity uq and the turbulent correlation length 
^0, is usually used in theoretical analysis. 

Figure [T^] demonstrates that a difference between the 
correlation functions determined in our experiments with 
and without transformation of images into a binary form 




FIG. 12: Normalized two-point second-order correlation func- 
tion $(-R) for aj, — 1 and different values of parameter a^: 
Uc = 1.2 (dashed-dotted line), = 2 (solid line) and Oc = 4 
(dashed line). Two-point second-order correlation function 
$(-R) determined in our experiments with (unfilled circles) 
and without (stars) transformation of images into a binary 
form. The distance R between two points is measured along 
the y axis ai z — (at the mean boundary of the jet). Here 
Re= 2 X 10^ and Pe= 2 x 10*. 
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FIG. 13: Two-point second-order correlation function ^{R) 
determined in our experiments without transformation of im- 
ages in a binary form. The distance R between two points is 
measured along the y axis: at the mean boundary {z — 0) of 
the jet (filled squares); at z = —0.55D inside the jet (filled 
triangles); at z — 0.55D outside the jet (stars). Here Re= 10* 
and Pe= 2 x 10*. 



is fairly small. Inspection of Fig. [T3l shows that the corre- 
lation functions are weakly dependent on z. This conclu- 
sion is also valid when we transform images into a binary 
form. Note that when we determine the two-point cor- 
relation function using images transformed into a binary 
form, we take into account the contributions of both, the 
jet fiuid and the ambient fluid. On the other hand, when 
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we determine the parameter A (see Sec. Ill), we take into 
account only the contribution of the jet fluid. The dif- 
ference between the results obtained using a measured 
random particle number density field and binarized (us- 
ing a threshold) random field is not surprising, since a 
binarization is equivalent to a nonlinear transformation 
(filtering) of the particle number density field. Note that 
the two procedures for data analysis of mixing, without 
and with transformation of images in a binary form, yield 
an additional information about mixing of particles at the 
external boundary of a jet. Solid line (a^ = 2) in Fig. 11 
approximates our experimental results after binarization 
reasonably well. 

This study demonstrates that there is a qualitative 
agreement between the measured and theoretically pre- 
dicted two-point second-order correlation functions (see 
Fig. [T^ . This is an indication that the presented theo- 
retical model of the particle number density fluctuations 
generated by tangling of the gradient of the mean par- 
ticle number density by the turbulent velocity field, can 
mimic mixing at the external boundary of a submerged 
turbulent jet. 

The compressibility of the turbulent diffusion tensor 
caused by the compressibility of Lagrangian trajectories 
in a turbulent flow with a finite correlation time is an 
universal feature of turbulent mixing. The direct conse- 
quence of the compressibility is an existence of a positive 
maximum in the potential U{R) in the scales R < Iq. 
This is a reason why the correlation function of the parti- 
cle number density does not have a power-law behaviour. 
This property of the correlation function of the particle 
number density is a generic feature, that seems to be 
independent of the origin of a random fluid flow. In par- 
ticular, a complicated random flow which arises at the 
nonlinear stage of evolution of Rayleigh- Taylor instabil- 
ity has a finite correlation time. This implies that the 
turbulent diffusion tensor in such flow is compressible. 
It is plausible to suggest that the latter is a reason that 
the PDF of the phase function found in the analysis of 
Rayleigh- Taylor instability in ^1^ does not have univer- 
sal scaling and cannot be approximated by a power law 
function. 



V. CONCLUSIONS 

In this study we investigate experimentally and theo- 
retically mixing at the external boundary of a submerged 
turbulent jet. An air jet in the experiments is seeded 
with the incense smoke particles having large Schmidt 
number and small Stokes number. The spatial distri- 
butions of the jet and ambient fluids are determined in 
these experiments using Particle Image Velocimetry and 
an Image Processing Technique based on the analysis of 
the intensity of the light Mie scattering. 

Two approaches are used in the data analysis, whereby 
one approach is based on the measured phase function for 
the study of the mixed state of two fluids. This approach 



is similar to that used for analysis of the Rayleigh- Taylor 
instability p^ . Other approach applied in this study 
is based on the analysis of the particle number density 
fluctuations generated by tangling of the gradient of the 
mean particle number density by the turbulent velocity 
field. This gradient is formed at the external boundary 
of a submerged turbulent jet. 

The main result of this study is that the phase- 
function-based formalism and the formalism based on the 
analysis of the two-point second-order correlation func- 
tion of the particle number density imply the absence 
of the universal scaling: the PDF of the phase function 
of a jet fluid penetrating into an external flow and the 
two-point second-order correlation function of the parti- 
cle number density exhibit a non-power-law behaviour. 

A power-law behavior of the correlation function of the 
particle number density can be possible when the char- 
acteristic scale of the inhomogeneity of the mean particle 
number density is larger than the maximum scale of tur- 
bulent motions. However, in the vicinity of the jet bound- 
ary the characteristic scale of the inhomogeneity of the 
mean particle number density is much smaller than the 
maximum scale of turbulent motions. It is plausible to 
suggest that one of the reasons for the non-universal be- 
havior of the correlation function of the particle number 
density is the compressibility of Lagrangian trajectories 
in a turbulent flow with a finite correlation time. For a 
delta-correlated in time random velocity field the power- 
law behavior of the correlation function of the particle 
number density might be possible. However, this sugges- 
tion should be carefully verified in different laboratory 
experiments and numerical simulations. 

A mechanism of mixing at the external boundary of a 
submerged turbulent jet is a kinematic tangling process. 
The tangling mechanism is universal and independent of 
the way of generation of turbulence for large Reynolds 
numbers. Tangling of the gradient of the large-scale 
velocity shear produces anisotropic velocity fluctuations 
40, 41] which are responsible for different phenomena: 
formation of large-scale coherent structures in a turbu- 
lent convection [i^], excitation of the large-scale inertial 
waves in a rotating inhomogeneous turbulence [1^, gen- 
eration of larg e-scale vorticity [11] and large-scale mag- 
netic field [43 in a sheared turbulence. 

Our theoretical study of mixing is based on the anal- 
ysis of the two-point second-order correlation function 
of the particle number density. There is a qualitative 
agreement between theoretical and experimental results 
obtained in this study. However, the theoretical results 
obtained in this study cannot be valid in the most general 
cases since we adopted a number of simplifying assump- 
tions about the turbulence. We considered a homoge- 
neous, isotropic, and incompressible background turbu- 
lence. We assumed also that the generated "tangling" 
fluctuations of particle number density do not affect the 
background turbulence. The latter assumption implies 
a one-way coupling between particles and fluid which is 
reasonable in view of a small mass-loading parameter for 
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particles. We adopted a simplified physical model in or- 
der to describe mixing of particles at the turbulent jet 
interface. In spite of the fact that the simplified model 
considered in our paper can only mimic mixing in real 
jets, this model describes robust features of turbulent 
mixing in our experiments. Clearly, the comprehensive 
theoretical and numerical studies are required for quanti- 
tative description of mixing at the turbulent jet interface. 
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APPENDIX A: ASYMPTOTIC SOLUTIONS FOR 
SECOND-ORDER CORRELATION FUNCTION 

In this appendix we perform an asymptotic analysis 
of solution for the two-point second-order correlation 
function ^(i?). The correlation function ^{R) satisfies 
the following boundary conditions: ^(i? = 0) = 1 and 
— > oo) = 0. This function has a global maximum 
at i? = and satisfies the conditions: ^'{R = 0) = 0, 
4>"(i? = 0) < and = 0) > > 0)|. There 

are several characteristic regions in the solution for the 
correlation function $(i?): (i) the viscous range < i? < 
Re"^/*; (u) the inertial range Re"^/"* < i? < 1 and (ui) 
the large scales R ^ 1 whereby there is no turbulence. 
Here R is measured in units of the maximum scale of 
turbulent motions ^o- 

The behaviour of the dimensionless two-point second- 
order correlation fimction $(-R) is described by Eq. pT|) . 
Equation (fTTj) with / = has an exact solution in the 
viscous range < i? < Re^"^^^: 



$(A) = Six) x-^ (1 + x^y/^ 



(Al) 



where S{X) is a real part of the function AiP^{iX) - 
A2Ql{iX). Here P^{iX) and Q|f(iX) are the Leg 
endre functions with imaginary argument, X^ = 
/3A/Pe-\/Rei?^, and 



C 

/3 



5ct 



1 + 3cr ' 
2 2 



1 + 3cr ' 



/5. 



M 



l + 3cr 
3(1 + cr^) 

4|7| 



/3mvRc 



1/2 



a = UcCT^ and 7 < 0. In order to obtain the solution (|A1|) 
we take into account that in the viscous range of scales 
the inverse mass, the function xi^) and the potential 
U{R) are given by the following formulas: 



U{R) = 



Pe 
2a- 



{1 + X' 



1-/3; 



x{R) = -^MiR)Re'/'R, 
X^ 



1 



1 + 



1 ^ 

/3i = 



hX2 
3(1- 



8a) 



l + 3a 



Here we use asymptotics of the functions F{R) and Fc{R) 
determined for R < Re"^/^ (see Eqs. ([241)- ([Ml) and pSjl - 
(dni)). The solution (|A1[) has the following asymptotics. 
In the scales < i? < l/[%/PeRe^/'*] (i.e., for X < 1), 
the solution for the correlation function $(i?) reads: 



$(i?) = 1 
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(A2) 



while in the scales l/[\/PeRe^/'*] < i? < Re"^/" (i.e., for 
X » 1), 



$(i?) = Ai + A2 R~ 



(A3) 



In the inertial range Re"^/* < i? < i?max, Eq. ^ 
with / = reads 



i?2$" + (l-/3)i?$' 



2M 

/3. 



i?2/3 $ = 



(A4) 



M 



where 



f3 



3 + 7a 



M 



2(3 7a) 

9(1 -hcr^) 



Here we use asymptotics of the functions F{R) and Fc{R) 
determined for Re"^/^ < i? < 1 (see Eqs. (|24l) -(|26 l) 
and (|12 p - (fT5|) ). Equation (|A4p has an exact solution: 



(A5) 



$(i?) = i?'^/2 j^(^) ^ Y,{z 



where J,^{z) and Y,y{z) are the Bessel functions of the first 
and the second kinds, z = P2 R^^^, v = (3/2) |/3| and /32 = 

— 1/2 

3(2|7|//3m) . Here we take into account that in the 
inertial range the inverse mass M-i(i?) = PmR^I'^ and 
the function x(-R) = -(4/3c//9)M(i?) i?!/^. The solution 
^5P has the following asymptotics in the scales Re^'^^'' < 



-1/3 



(i.e., for z <^ 1): 

^{R) = As + AiR-^ 



(A6) 



In the vicinity of the maximum of the potential U 



Umax that is attained at the scale Rr, 



1, the inverse 



mass ~ 1 and it changes with R considerably slower 
than in the range R <ti 1. The asymptotic solution for 
the correlation function <&(-/?) in this region reads: 



HR) = As 



CXp [ - k{R - i?max)] 

RP 



(A7) 
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where k = (C/max — 2|7|) and & ^ 1. 

For large scales, 1, there is no turbulence so that 
~ 2/3 and x = U = 0. The asymptotic solution 
for the correlation function $(-R) in this range of scales 
reads: 

= AqR-^ sin [R ^3171/2 + tpo) • (A8) 
Matching functions $(-R) and ^'{R) at the boundaries 



of the above-mentioned regions yields coefficients Ak and 
the damping rate 7. In particular, Ai ~ A3 ~ ^5 ~ 
Aq fa 1, and A2 ^ 1, A4 ^ 1, where we consider the 
case Sc^ 1. The results of the above asymptotic analy- 
sis are in agreement with the numerical solution for the 
correlation function $(-R) obtained in Sec. IV. 
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